Ticking terahertz wave generation in attoseconds 
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We perform a joint measurement of terahertz waves and high-order harmonics generated from 
noble atoms driven by a fundamental laser pulse and its second harmonic. By correlating their 
dependence on the phase-delay of the two pulses, we determine the generation of THz waves in tens 
of attoseconds precision. Compared with simulations and models, we find that the laser-assisted 
soft-collision of the electron wave packet with the atomic core plays a key role. It is demonstrated 
that the rescattering process, being indispensable in HHG processes, dominant THz wave generation 
as well but in a more elaborate way. The new finding might be helpful for the full characterization 
of the rescattering dynamics. 



Intense terahertz (THz) radiation from ionizing gases 
in two-color laser fields has attracted great interest re- 
cently not only for the application such as remote sensing 
[l[ but also for understanding the mechanism of terahertz 
generation 0413 ■ When the fundamental and frequency- 
doubled femtosecond laser pulses are focused in gases, the 
intensity of the THz radiation is modulated by the rela- 
tive phase between the two pulses. Experiments confirm 
that THz output increases dramatically with the onset of 
plasma formation, indicating that ionization plays a key 
role in THz generation 0, E[ • This led to involve four- 
wave mixing (FWM) in plasma to elucidate the extreme 
nonlinear phenomenon but without the knowledge of the 
origin of v^ 3 ' responsible for such high efficiency. 

A transient photoelectron current model is recently 
proposed taking into account of the electron dynam- 
ics following ionization [H, 0]- In this model, electrons 
are released from atoms through tunnel ionization and 
then driven by non-sinusoidal two-color laser field which 
breaks the dynamic symmetry of the electronic wave- 
packet. A directional current is thus formed causing the 
emission of THz waves. However, there is a discrepancy 
between the above models on the relative phase when 
the terahertz yield maximizes. It is evident that the 
key problem of describing THz generation in two-color 
fields is to find the buildup process and the subsequent 
relaxation of the transient current 0, [H| ■ The sensitiv- 
ity on the phase delay indicates that the sub-cycle elec- 
tron dynamics needs to be measured to understand THz 
generation and a determination of the relative phase in 
attosecond precision is urgently expected. 

It is known that a commensurate two-color field can be 
used to steer trajectories of continuum electrons and con- 
trol their re-scattering/re-collision with the parent ion to 
manipulate high-order harmonic generation (HHG) f]~3l — 
Il5j |. In HHG processes, each harmonic can be uniquely 
related to two quantum orbits within each optical cy- 
cle based on the rescattering theory [l6|. Two-photon 
ionization experiments have confirmed that the emission 
time of harmonics differs in attoseconds and it varies lin- 
early with harmonic order in positive slope for the short 
trajectories 17-13]. The weak harmonic field tunes the 



phase difference of radiation in the positive and negative 
half cycle of the fundamental, causing the appearance of 
even-order harmonics with intensity modulated by phase 
delay of the two pulses. The two-color scheme thus pro- 
vides an all-optical in situ method for measuring the birth 
of XUV pulses [l3j which is an ideal tool for ticking the 
generation of THz in attoseconds. Furthermore, since 
both yields of THz waves and harmonics modulate when 
the time-delay of the two pulses varies by attoseconds, a 
question arises: are they related knowing that THz waves 
have duration of picoseconds? 

We therefore perform a joint measurement by monitor- 
ing THz and high-order harmonic yields in a two-color 
laser field simultaneously. By correlating their phase- 
delay dependence, the absolute relative phase depen- 
dence of THz yields can be retrieved and it is found 
that THz yields take maximum at phase delay of 0.87T 
deviated from the prediction of previous models. Com- 
pared with numerical simulations and theoretical models, 
we find that the Coulomb potential from the residue ion 
should be taken into account and the laser-assisted soft- 
collision of the electron wave packet with the atomic core 
affects the phase dependence of Thz yields. Thus tera- 
hertz radiation provides an complementary tool to study 
the electron dynamics in addition to high-order harmonic 
generation and above threshold ionization, helping the 
full characterization of the rescattering wave packet. 

In the experiment, we used a 1-kHz ThSapphire 
chirped pulse amplification (CPA) system that delivered 
1.6 mJ at 790 nm in a 25 fs pulse. The beam is re- 
duced to a diameter of 10 mm and split into pump and 
probe beams. The pump is focused on a 200 /im thick 
/3-barium borate (BBO) crystal (type I) to generate the 
second harmonic beam with 20% conversion efficiency. 
A two-color Mach-Zchndcr interferometer is employed to 
control the relative phase between the two pulses, with 
a pair of fused-silica wedges in the fundamental path. 
Neutral density filters are used in both arms to control 
their intensities separately. High-order harmonics and 
THz waves are generated in the 1.5 mm long gas cell in 
12 mbar by the two-color pulse using an off-axis parabolic 
Au-coated mirror with focal length of 200 mm. The gas 
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cell is placed 2.5 mm behind the focus to optimize the 
HHG yields. The THz waves are reflected with a hole- 
drilled off-axis parabolic mirror and then combined with 
the delayed probe beam into a 1mm thick (110) oriented 
ZnTe crystal for electro-optic sampling to the waveform 
below 3 THz. The harmonics passing the hole is simulta- 
neously recorded by a home-made flat-field spectrometer 
for each phase-delay. 

In Fig. [1] we show the measured THz wave form and 
the yields of high-order harmonics modulated with the 
time delay of the two-color laser field which is charac- 
terized by E(t) = Ei cos ujt + E2Cos(2uit + </>), where 
the field strengths of the fundamental and the second- 
harmonic are denoted by E\ and E2 respectively, and 
the phase delay is related to the time delay by <f> = 2cutd- 
The solid line depicts the phase-delay tf>o(Nu) maximiz- 
ing the yields of even-order harmonics from 16th to 26th. 
In order to determine the absolute phase-delay, we apply 
the procedure given in jl3} by mapping <po(Nui) to the 
emission time of harmonics obtained from the rescatter- 
ing model [l6 . The intensity of 2N th order harmonic is 
modulated approximately as sin 2 A, where A is the ad- 
ditional action introduced by the vector potential of the 
second-harmonic pulse and is given by [1 31 ) 

Ass — \ I (sin cut — sin lotA sin(2wi + <j))dt (1) 
2uj hi 

with the ionization instant ti and the recombination 
instant t r as functions of harmonic energy and A is 
the ratio of E 2 to E\. The corresponding phase-delay 
that maximizes 2N" 1 harmonic yields is plotted in solid 
lines in Fig. [2j The laser intensity is estimated of 
2.3xl0 14 W/cm 2 from the cut-off energy of harmonics. 
Also shown in Fig. [5] is the extracted experimental data 
denoted in circles. The lower order harmonic is tak- 
ing maximum at later delay in accordance to the har- 
monic generation from short trajectories Note that 
the left part of the theoretical result (solid line) corre- 
sponds to the case of harmonic generation from the long 
trajectories which has been removed by phase-matching 
in experiments. The overall agreement between theory 
and experiment for all the harmonics presented allows us 
to determine the absolute phase-delay maximizing yields 
of harmonics, e.g., the 22rd harmonic takes maximum 
at 0.937T. Note less than 0.047T phase-delay control is 
achieved in the experiment corresponding to 60 attosec- 
onds. 

Now using the harmonic modulation as ticks we can as- 
sign the absolute phase-delay to the modulation of THz 
yields (summed yields below 3 THz) shown in Fig. [2] Be- 
cause the THz and harmonic yields are measured jointly 
for each given phase-delay, we can in fact lock the mod- 
ulation of THz yield to one of the harmonics. The THz 
yields are found taking maximum at phase delay of 0.87T. 
Using the quantum mechanical approach 0, [2Cfl , we ob- 
tain the same result which deviates from the predication 
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FIG. 1: Contours of the terahertz wave forms (bottom) and 
the intensities of harmonics (top) from 16th to 26th varying 
with the time delay of the two-color laser field, in unit of 
the fundamental optical period. The intensity ratio of the 
second-harmonic to the fundamental field is about 0.5%. 




Phase-delay (n) 

FIG. 2: Top: The phase-delay maximizing intensity of even- 
order harmonics extracted from the experiment (circles) and 
theoretical model (solid line). Bottom: modulation of THz 
yields with phase-delay. The laser intensity is estimated of 
2.3xl0 14 W/cm 2 . 

of four-wave-mixing and directional photoelectron cur- 
rent models. What is the reason behind? 

According to the photoelectron current model, the mi- 
cro current (j(U)) generated from the ionization event at 
time ti can be obtained by weighting the final outgoing 
electron velocity with the corresponding instantaneous 
ionization rate @. The initial velocity of the electron 
born at time ti is assumed zero and the final electron 
velocity can be found from the classical motion to the 
electron in laser fields. The total current responsible for 
THz emission is thus the coherent summation over all the 
ionization events. Because the micro current generated 
from two ionization events separated by tt/lo has symme- 
try of j(ti) — —j(U +tt/(jo), their destructive interference 
causes vanishing total currents for THz generation in a 
single color laser field. When a second harmonic pulse 



3 



is introduced, the symmetry between tj and tj + tt/lo is 
broken, directional current is formed giving rise to the 
emission of THz waves. 

Depending on the instants it is set free, the electron can 
either directly escape from the atom or recollide with the 
atomic core. We can therefore define the escaping cur- 
rents (EC) and rescattering currents (RC), j eS c{U) and 
jres(ti) accordingly. Because of the large energy gained 
from the laser field by electrons approaching the core, 
HHG is mainly contributed by the high frequency com- 
ponents of the RC. On the contrary, THz wave has pho- 
ton energy in meVs, it can be emitted at much large 
distances from the nucleus involving free-free transitions 
[20j|. Therefore the contributions from low energy and 
less accelerated electrons can not be ruled out and the 
EC needs to be considered as well. Note that EC/RC 
is related to tunnel ionization occurring in the quarter 
cycles before/after the electric field maxima. Because 
jresiU) — —jesc(—ti), the magnitude of the total cur- 
rent from a two-color laser pulse takes maximum at the 
relative phase of n/2 0]- 

In the above discussion regarding free electron mo- 
tion, only the electron-light interaction is considered and 
the influence of the Coulomb potential from the residue 
atomic core is neglected completely. This so-called strong 
field approximation (lif has been widely used in modeling 
processes such as HHG where high kinetic energy elec- 
trons are involved. However for near-threshold phenom- 
ena, it can distort the free-electron wave function signifi- 
cantly and modify the density of states 2l| . Its focusing 
effects enhance double ionization by one order of magni- 
tude [2^ |. The Coulomb- laser coupling has been shown 
affecting laser-assisted photoionization ellipticity- 
dependence of low order harmonics [j^ j. angular distri- 
bution [25| and low-energy peaks of of above-threshold 
ionization (ATI) photoelectrons 26[ as well. Since THz 
wave generation involves of slow electron motion which 
can be strongly distorted by the Coulomb potential, the 
response of the whole wave packet in the combination of 
the laser field and the Coulomb potential needs to be in- 
vestigated. In the following we show that the Coulomb 
potential has different impacts on the formation of the 
EC and RC. 

We consider two ionization instants t\ and ti in the 
two-color laser pulse with zero phase-delay. The cor- 
responding classical trajectories of the electron driven 
by the laser field only are shown in the dotted lines 
in Fig. |3] (a) and (d) respectively. Quantum mechani- 
cally, we propagate the continuum electron wave packet 
(EWP) following tunneling ionization. The initial elec- 
tron wave packets at the instants of tunneling ionization 
are assumed to have the same shape of the ground state 
wave function but shifted to the classical turning points 
at xq = —Ip/Ef, where I p is the ionization potential 
and Ef is the corresponding instantaneous electric field. 
One-dimensional time-dependent Schrodinger equation is 





FIG. 3: Time propagation of electron wave packets after tun- 
neling ionization at two instants t\ (left column) and ti (right 
column) corresponding to the cases of escaping and rescat- 
tering respectively. The contours plots show the normalized 
densities of probability as functions of time and coordinates, 
encoded with color scale 1.0 corresponding to the density of 
0.02 a.u.. In (b) and (e), the Coulomb potential is neglected, 
but being considered in (c) and (f). In (a) and (d), the dotted 
lines represent the classical trajectories of electrons with zero 
initial velocity, and the solid lines represent the quantum me- 
chanical averaged electron displacements. The intensities of 
the fundamental pulse and its second harmonic are 1 x 10 14 
W/cm 2 and 5 x 10 12 W/cm 2 respectively, with zero phase de- 
lay. Both pulses have durations of 25 fs with Gaussian envelop 
centered at zero time. 



solved for simplicity. The evolving of the EWP's in the 
absence of Coulomb potential are shown in figure. G^b) 
and (e) respectively. As one of the quantum features, 
the EWP's are broadened during evolution because of 
the quantum dispersion. However, the center motion of 
the EWP coincides exactly with the classical motion of 
the electron in a laser field with zero initial momentum, 
thanks to the conservation of the canonical momentum 
in both the classical and quantum theories. 

It can be seen in Fig. [3Je) that the wave packet moves 
toward the negative X axis initially and then turns back 
to the positive axis passing the origin where the atomic 
core is positioned. It is clearly the scenario of rescat- 
tering. The central part of the wave packet makes hard 
collision with the atomic core that causes high-oder har- 
monic generation and high energy above-threshold ion- 
ization when followed by either recombination or back- 
scattering [27}. However, the THz yields are determined 
by the ejecting photoelectron current, therefore we must 
consider the propagation of the wavepackt after scatter- 
ing by the atomic potential. We see in Fig. [3] (f) the 
electron wave packet is fragmented after re-collided with 
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the atomic core. Especially the upper wing of the elec- 
tron wave packet in Fig. [3] (f) is strongly enhanced in 
density comparable to that of the center of the electron 
wave packet. The magnitude of the averaged velocity is 
increased because of more fractions of electron moving 
toward the positive direction. The averaged displace- 
ment plotted in solid line in Fig. [3jd) demonstrates the 
photoclcctron current is indeed enhanced. Different from 
recombination and backscattering processes, here we are 
considering the evolution of the whole wave packet with 
most of the packet scattered away from the atom. In 
other words, it is the laser-driven soft collision of the 
EWPs with atomic core causing the enhancing of the 
rescattering current. It might be the same mechanism 
that causes low-energy peaks in photoelectron spectra 
& 

On the contrary in Fig. |3fb) , the center of the EWP 
moves away from the origin in a manner of quiver mo- 
tion, without the chance to reencounter the atomic core 
making hard collisions. Due to the broadening of the 
EWP, fractions of electrons in the wing of the EWP pass 
through the origin, but the averaged velocity of the wave 
packet makes no change when the atomic potential is ig- 
nored. If the Coulomb potential is included in the time 
propagation we see in Fig. [3tc) that the wing of the EWP 
is altered by the atomic core. The EWP is fragmented as 
well with more fractions of electrons moving toward the 
positive axis by the attractive Coulomb potential. Be- 
cause the total drift velocity is negative in the absence 
of Coulomb potential, it is suppressed as plotted in the 
solid line compared with the dotted line in Fig. [3] (a). 
Here we see again that the soft-collision with the atomic 
core has profound effect on the EC. Note EC makes no 
contribution in the theory of HHG [n| . 

The laser-assisted soft collision is mainly occurring 
within the first quarter cycle following ionization because 
the electron drift away from the nucleus at later time such 
that the effect of Coulomb potential is not significant any 
longer. The momentum transfer is contributed most from 
the time that the electron makes the turn where its ve- 
locity is zero such that it has more time interacting with 
the Coulomb potential. It can be investigated theoreti- 
cally from the classical electron motion in the combina- 
tion of the atomic potential and the laser field 23 , 29| . In 
the scenario of rescattering, our estimation gives approx- 
imate momentum transfer of Ap 2E cos 8/ k 3 , where 
k = \J21 V with I p of the atomic ionization potential and 
E cos 8 is the field at instant of ionization. Because of the 
momenta transfer by the Coulomb potential, the sym- 
metry between the EC and the RC is broken with the 
former suppressed and the RC enhanced. Therefore the 
rescattering process responding for HHG, dominates THz 
wave generation as well but because of the laser-driven 
soft-collision of EWP with atomic core. 

In order to further check the analysis above, we per- 
form classical and quantum calculation for the THz yields 




FIG. 4: (a) Escaping, rescattering and total currents calcu- 
lated from the photocurrent model and (b) the corresponding 
THz yield below 3 THz; (c) THz yields obtained from classi- 
cal wavepacket simulation by considering Coulomb potential; 
(d) THz yields calculated from solving the time-dependent 
Schrodinger equation with the same laser parameters given in 
Fig. O 



from a two-color laser pulse. As shown in Fig. Hta) and 
(b), the THz yield takes maximum at phase-delay of 0.57T 
from the photocurrent model which ignores the Coulomb 
potential. When the Coulomb potential is taken into 
account in the classical wave packet simulation shown 
in Fig. |Uc), the THz yield from the EC is suppressed 
and that from the RC is enhanced. In particularly, the 
total THz yield takes maximum at the phase-delay of 
0.87T agreed with one dimensional quantum mechanical 
calculation shown in Fig. [4] (d) confirming the previous 
discussion. In the classical wave packet simulation, we 
sample and trace the trajectories starting with Gaussian 
distributed transversal momenta according to the tunnel- 
ing ionization theory [3(ij |. 

Finally, the phase-delay dependence of the RC can be 
obtained analytically by treating the second harmonic as 
a perturbation. It yields (details will be given elsewhere) 
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where e and m represent the charge and the mass of the 
electron respectively, and a — 2k 3 /3_Ei, where ionization 
e~ Q for hydrogen-like atom is used [30]. The 
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generated THz wave field strength is oscillating with <j) 
at the period of 2-7T which corresponding to time delay 
of 7r/w. It can be seen that THz wave takes maximum 

when tan <j) = ■ This critical phase is mainly de- 

termined by the ionization potential and the intensity of 
the fundamental laser pulse. For weak E\ when a> 1, 
the maximum yield occurs at zero phase agreed with the 
perturbative FWR theory. For the intensity of the fun- 
damental we uses, the maximum occurs at phase of 0.87T 
agrees with our measurements. For intensity much higher 
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the phase approaches to 7r/2 in accordance to the tran- 
sient current model [7J. Therefore the two mechanisms 
are channelled depending on the competition of the soft- 
collision and the electron-light scattering. 

In conclusion, the joint measurement of THz and har- 
monic yields allows us to determine the generation of THz 
waves in attosecond precision. While the time scales of 
the two kinds of radiation are quite different, we show 
that in fact they are connected and both originate from 
the same electron dynamics: laser-driven electron mo- 
tion in the atomic fields. Particularly we found that THz 
generation is sensitive to the Coulomb potential of the 
residue ion after tunneling ionization. The magnitude 
of current from directly escaping electrons is suppressed 
while that from the rescattering electrons is enhanced due 
to the soft-collision of electron wave packet with the ion. 
As a first step toward fully understanding the correlation 
of THz generation and HHG, we show that rescattering 
electron dynamics might bridge the two radiation pro- 
cesses. While HHG reflects the motion of the center of 
electron wave packt which makes hard collision with the 
atomic core, THz generation reflects the whole dynamics 
of the electron wave packet. Therefore the THz wave pro- 
vides a tool to investigate the whole electron wavepackt 
dynamics, beyond the recombination in HHG and back- 
scattering in above threshold ionization. We expect that 
the sensitivity of THz yields on Coulomb potential could 
be used to map atomic fields from within and help the 
full characterization of the rescattering wave packet. 
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